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Abstract. Galactic disks can form in asymmetric potentials of the assembling dark matter 
(DM) halos, giving rise to the first generation of gas-rich bars. Properties of these bars differ 
from canonical bars analyzed so far. Moreover, rapid disk growth is associated with the 
influx of clumpy DM and baryons along the large-scale filaments. Subsequent interactions 
between this substructure and the disk can trigger generations of bars, which can explain 
their ubiquity in the Universe. I provide a brief summary of such bar properties and argue 
that they fit naturally within the broad cosmological context of a hierarchical buildup of 
structure in the universe. 
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1. Introduction 

Among various approaches to galaxy forma- 
tion, two stand out — those that deal with iso- 
lated and open systems. In the former modus 
operandi, angular momentum, /, and mass of 
a system are conserved. Thus, for example, one 
deals with a finite reservoir of gas, galaxy inter- 
actions are treated in maximally-controled ex- 
periments, and, in most such cases, the amount 
of substructure present is quite limited, if it ex- 
ists at all. The modeled galaxies quickly use up 
the available gas and show a single splash of 
gas influx to fuel the central starburst or AGN 
activity. As a result, there are numerous diffi- 
culties in reproducing the observed phenomena 
in isolated systems: they start with simplified 
initial conditions, e.g., an axial symmetry (in 
disk galaxies) which is broken spontaneously 
by the bar instability. In many ways these sys- 
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terns are too idealized to explain the evolution- 
ary tracks of real galaxies, but they serve as ex- 
cellent educational tools. The perennial ques- 
tion about these systems remains: why should 
nature create 'equilibrium' systems and imme- 
diately 'destroy' them by various dynamical 
instabilities? 

Here we shall take an alternative approach 
and treat an assembling disk-halo as an open 
system whose evolution is driven by both ex- 
trinsic (i.e., environmental) and intrinsic fac- 
tors. In this context, 'open' is synonymous with 
'cosmological.' The cosmological caveats in- 
clude the redshift dependence of gas and dark 
matter (DM) influx along the large-scale fila- 
ments which feed the growing disk inside an 
assembling halo. Thus, conditions which are 
associated with a disk-halo growth are those of 
an asymmetric background DM potential and 
clumpy mix of DM and baryons. This asymme- 
try of the potential comes from the surrounding 
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Fig. 1. Snapshots of gas density (color-coded) evolution in the BDM model at z = 9 (based on Romano- 
Diaz et al. 2008a). The frames are from left-to-right 600 kpc, 200 kpc, 50 kpc and 25 kpc across. 



filaments(s) and from virialized triaxial DM 
mass distribution. The accretion process itself 
is highly inhomogeneous and potentially re- 
sponsible for various dynamical phenomena in 
the disk and its immediate surroundings — few 
of these will be discussed here. 

Galactic bars appear as dominant morpho- 
logical features in nearby disks (e.g., Sellwood 
& Wilkinson 1993; Knapen et al. 2000; 
Marino va & Jogee 2007). At intermediate red- 
shifts of z ~ 0.2-1, this picture is a subject 
to controversy (Jogee et al. 2004; Elmegreen 
et al. 2004; but see Sheth et al. 2008). At even 
higher z we know very little about the disk evo- 
lution. However, the importance of bars goes 
well beyond morphology — they serve as a 
fundamental channel of angular momentum re- 
distribution in a galaxy. Hence bar formation 
and evolution should be tied to the overall pro- 
cess of galaxy formation and evolution in the 
CDM universe. Here we attempt to pursue this 
argument. 

2. First Generation of Bars in Open 
Systems 

Cold stream-fed growth emerges as the dom- 
inant mode of DM and baryonic assembly in 
massive galaxies (e.g., Dekel et al. 2008). In 
the following discussion, we rely on represen- 
tative models which run from identical initial 
conditions at z = 120 in ACDM Universe — 
pure DM (PDM) simulations and simulations 
with WMAP3 fraction of baryons, the BDM 
models (Romano-Diaz et al. 2008a). Fig. 1 dis- 
plays snapshots of gas density evolution in the 
BDM model on scales from 600 kpc to 25 kpc. 



The latter scale is relevant for disk formation 
and it is significant that the filamentary struc- 
ture extends down to this size at high z. The 
resulting halo shapes will be substantially tri- 
axial (e.g., Allgood et al. 2006), partly because 
of the assembly process and partly as a result 
of the radial orbit instability. Baryons act as 
to wash out the equatorial ellipticity and to re- 
duce the flatness of the DM distribution (Fig. 2; 
also e.g., Kazantzidis et al. 2004; Berentzen & 
Shlosman 2006; Shlosman 2007). 

One of the crucial differences between the 
local and high-z universes lies in that the lat- 
ter is much more cold gas-rich. Gas-dominated 
disks that form in the asymmetric potentials 
cannot maintain axial symmetry — stars and 
especially gas will respond dramatically pro- 
ducing a standing density wave and a shock, re- 
spectively. Simulations have demonstrated that 
early bars will form under these conditions 
and that these bars will have properties dif- 
ferent from the canonical bars discussed in 
the literature so far (e.g., Heller et al. 2007b). 
Specifically, this first generation of bars (1) 
will form in disks with a low surface density 
and central mass concentration. Consequently, 
(2) they will exhibit low pattern speeds. (3) 
They will experience a rapid growth, concur- 
rently with that of the host galactic disks. 
Finally, (4) both bars and disks will be aware 
of the DM halo orientation, i.e., position of its 
major axis. Importantly, the halo figure tum- 
bling is exceedingly slow, whether it grows in 
a filament or in the field (Bailin & Steinmetz 
2004; Heller et al. 2007b; Romano-Diaz et 
al. 2008a,c). This will result necessarily in a 
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Fig. 2. Evolution of the axial ratios (shapes) for the prime DM halo at various representative radii in the 
PDM (left) and BDM (right) simulations. Top: halo b/a (stratified roughly from top to bottom) at r = 10 kpc, 
20 kpc, the NFW radius, radius of maximal circular velocity, 100 kpc, 200 kpc, 300 kpc and the virial radius. 
Bottom: the same for c/a (Romano-Diaz et al. 2008a). 



double inner Lindblad resonance (ILR) — the 
outer ILR will move far out and the inner ILR 
to the center — between these resonances the 
baryonic response in the disk will be orthogo- 
nal to the major axis of the halo, reducing its 
equatorial ellipticity. 

Thus, first bars are expected to form in the 
growing disks embedded in the triaxial halos. 
However, bars and triaxial halos are incompat- 
ible — a canonical bar tumbling within a stag- 
nating halo will serve as a source of chaos and 
dissolve quickly (El-Zant & Shlosman 2002). 
The caveat lies in that the early bars are sup- 
posed to rotate very slowly (see above) and 
so the amount of chaos will be greatly re- 
duced, prolonging the bar dissolution time to 
-1-2 Gyr. During this time period, the disk 
growth can be appreciable. If the disk response 
to the halo forcing (see above) can success- 
fully wash out the inner halo equatorial ellip- 
ticity over this timescale, it will remove the 
main source of chaos and open the possibility 
for the bar to strengthen. Future work can con- 
firm whether this indeed what is happening. 

Bars living in growing disks are necessarily 
gas-rich (e.g., Heller et al. 2007 a,b; Romano- 
Diaz et al. 2008a,c), to the extent that was not 
contemplated in the literature so far. Shlosman 
& Noguchi (1993) found that, for gas fractions 
/ g > 10% of the disk mass, the Jeans insta- 



bility will lead to massive clumps which will 
spiral down to the central kpc as a result of a 
dynamical friction, leading to a nuclear star- 
burst. Heating the stellar disk on their way, the 
clumps can damp the canonical bar instabil- 
ity or postpone it substantially. Bournaud et al. 
(2005) claimed that f g ~ 7.25% is sufficient to 
trigger the bar dissolution due to the / trans- 
fer from gas to stars in the disk, in simulations 
with a rigid halo. Berentzen et al. (2007) com- 
pared the evolution of disks with f g ~ 0%-8%, 
concluding that the apparent decrease in the 
bar strength is not related to the gas content 
but rather to the vertical buckling instability 
(e.g., Combes et al. 1990; Raha et al. 1991) in 
the bar, which weakens it substantially but al- 
lows for a secular rebuilding (e.g., Martinez- 
Valpuesta et al. 2004; 2006). Bars that are 
capable of growing typically appear as 'fast' 
bars, i.e., extending to near corotation, even in 
cuspy halos (Dubinski et al. 2008). 

Gas fractions to be encountered at high 
Z can be an order of magnitude higher than 
quoted above. Results, based on the equilib- 
rium models, show that dependence of the bar 
strength on f g is not so dramatic. As long as 
the Jeans instability in the gas is damped, / 
transfer from the gas to the stars is compen- 
sated by / lost by the stellar bar to the DM 
halo. Fig. 3 demonstrates that development of 
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Fig. 3. Evolution of the stellar bar strength, given by the Fourier amplitude A 2 of the m = 2 mode, for 
various gas fractions, f g ~ 0% - 50% in the disk immersed in the live halo. The time is given in Gyrs. 
The Jeans instability is progressively damped from top to bottom models. The first peak defines the time 
of the vertical buckling in the bar. Note that A 2 shapes are very similar in all models up to the first peak 
(Villa- Vargas et al. 2008). 



the canonical bar instability is nearly indepen- 
dent of the gas fraction, up to f g ~ 50% at least, 
if the Jeans instability in the gas is suppressed 
(Villa- Vargas et al. 2008). The buckling pro- 
ceeds similarly independent of f g . But the sub- 
sequent evolution of the bar strength differs 
— while 15% gas still allows for the secular 
growth of the bar, 30% reduces the bar strength 
by a factor of 2-3 over the Hubble time. 

3. Dynamics of Bars in Open Systems 

When gas rich bars are triggered by the asym- 
metric DM potential, they can appear simulta- 
neously on more than one spatial scale. By this 



we mean that bars on different scales will either 
have different pattern speeds or will be situated 
at different position angles (in fact, orthogo- 
nal to each other) and have the same speed. 
Most frequently they form on smaller scales of 
a fewx 100 pc — so-called nuclear bars, and on 
the scales of a few kpc, so-called prime bars, 
although in some cases only one bar forms 
(Heller et al et al. 2007 a,b). As we have stated 
above, these bars do not form as a result of the 
canonical bar instability but rather constitute 
the disk response to the finite external pertur- 
bation. Unlike most of the bars in the local uni- 
verse, their young counterparts appear to be de- 
lineated by star formation which is present not 
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only close to the ILR(s) but throughout the area 
of the bar. In a way they resemble inclined star- 
bursting disks, unless distinguished kinemati- 
cally. 

As a first step, we discuss the behavior of 
single bars formed under these conditions, then 
comment on the evolution of double bars. Bar 
dynamics can be studied by inspecting their 
pattern speed, Db, evolution. Here again the 
'cosmological' bars differ from their low z gas- 
poor representatives. Those that are triggered 
by the asymmetric potential are characterized 
by strongly growing Db, concurrently with the 
disk growth. An alternative way of triggering 
bars, also by a finite perturbation, was known 
for some time as a tidal triggering (e.g., Byrd 
et al. 1986; Noguchi 1987), but was never con- 
sidered to be important enough to explain the 
observed bar fraction in the universe. As we 
discuss below and in the next section, within 
a somewhat modified framework, this picture 
may become more appealing. 

It is important to emphasize here that the 
class of tidally triggered bars exhibit initial 
Qb which are determined by the orbital speeds 
and the impact parameters of the perturbers. 
But the subsequent evolution of Qb in all three 
types of bars (namely, canonical, triggered by 
the DM asymmetry or triggered tidally) is de- 
termined by the balance of angular momentum 
in the bar and by the central mass concentra- 
tion, i.e., the shape of the inner rotation curve 
in a galaxy. Here, the role of gas, the only dis- 
sipative agent, cannot be overestimated. Either 
dynamically or secularly, the gas influx to the 
center contributes to the mass growth there and 
so, even under conservation of / in the bar, will 
cause an increase in Db- The net effect on the 
bar tumbling can be estimated from the amount 
of / supplied by the disk gas to the bar within 
its corotation radius, the growth of the central 
mass concentration, and the amount of / lost 
by the bar to the outer disk and especially to 
the DM halo. 

Simulations of isolated disks with an ad 
hoc addition of gas (Bournaud & Combes 
2002) and cosmological simulations of grow- 
ing disks embedded in large-scale filaments 
(Romano-Diaz et al. 2008c) demonstrated that 
bars can accelerate over a prolonged time pe- 



riod of a few Gyrs, during which the sur- 
face density of the disk, especially of its gas 
component, grows. This underlines the im- 
portance of the cold gas influx via filaments. 
With this, the strength, A 2 , of the cosmologi- 
cal bars (anti)correlates with Qb — weakening 
bars accelerate, and vice versa, as shown first 
for the pure stellar bars in equilibrium models 
(Athanassoula 2003). Exceptions include the 
time periods of interactions between the devel- 
oped bars and flybys, and the double bar sys- 
tems (Romano-Diaz et al. 2008c). 

In the case of double bars, only the nuclear 
bars accelerate after formation, while the pri- 
mary bar pattern speeds are quite steady over 
long time periods (Heller et al. 2007a). Both 
shapes and pattern speeds of nuclear bars vary 
with their orientation to the prime bar (Heller et 
al. 2001; Shlosman & Heller 2002; Englmaier 
& Shlosman 2004; Juntai & Debattista 2007) 
and this is visible in models of cosmological 
double bars as well (Heller et al. 2007a). The 
latter work has also shown that double bars can 
lock their pattern speeds at various ratios of 
nuclear-to-primary Q n /^p — an issue of mode 
coupling. 

4. Bar Triggering by DM Substructure 

The possibility of bar triggering via interac- 
tions with other galaxies hints at a strong de- 
pendence of bar fractions on the environment 
and redshift. However, bar fractions in clusters 
of galaxies and in the field appear comparable 
(Marino va et al. 2008; and these Proceedings). 
This discrepancy between expectations and ob- 
servations can be erased if bars can be trig- 
gered by interactions with DM substructure 
and not by a more rare interactions with their 
more massive neighbors. Existence of such 
DM subhalos is inherent to the hierarchical 
scenario of galaxy formation (e.g., White & 
Rees 1978). Hence, whether a halo is found in 
the cluster or in the field, it will grow by con- 
suming the surrounding subhalos — those that 
survive the radial plunge can have a profound 
effect on the disk evolution. 

Properties of a subhalo population have 
been studied in pure DM simulations (e.g., 
Ghigna et al. 1998; Klypin et al. 1999; 
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Fig. 4. Mass function A/^(M sbh ) of the subhalo population at z = 2 in PDM (blue dashed) and BDM (red 
dotted) models. The minimal resolved subhalo mass is 10 8 M Q , which is about 10~ 4 3 of the prime halo virial 
mass at z = 0. Left frame: N(M shh ) for subhalos outside the virial radius of the prime halo; Right frame: 
N(M shh ) for subhalos inside the central ~ 100 kpc x 50 kpc virialized region (defined by the isodensity 
contours) of the prime (Romano-Diaz et al. 2008d). 





Fig. 5. Population of subhalos in PDM and BDM models within the inner 60 kpc (left) and inner 30 kpc 
(right) of the prime halo at lower redshifts (Romano-Diaz et al. 2008b). 



Kravtsov et al. 2004; Diemand et al. 2007). 
Their mass distribution function was found to 
be nicely approximated by a power law Af S bh ~ 
M~£ h , where M s bh is the mass of a subhalo (e.g., 
Diemand et al. 2004). Much less is known 
about subhalos in the presence of baryons. 
Their power law distribution seems to be re- 
tained for DM halos above 10 10 M© — the res- 
olution limit (Weinberg et al. 2008). Fig. 4 dis- 
plays the subhalo mass functions for PDM and 
BDM models at z ~ 2, for much higher res- 
olution simulations resolving subhalos above 
10 8 M Q (Romano-Diaz et a;. 2008d). To em- 
phasize the environmental effects, the popula- 
tion was divided into the outer one, situated in 
the computational box outside the virial radius 



of the prime halo, and the inner one, situated 
within the central virialized region delineated 
by density isocontours, -100 kpc x 50 kpc. A 
number of trends can be seen: (1) for the 'field' 
subhalos, the difference between the PDM and 
BDM models is minimal at these redshifts. (2) 
the log-log slope is about (-0.9) - (-0.95) and 
quite stable at other z. (3) Evolutionary trends 
are present in the inner subhalo population — 
their integrated number is lower by a factor of a 
few, both in the PDM and BDM. (4) the BDM 
slope is shallower than the PDM one — more 
massive BDM subhalos survive better than the 
PDM ones in the central dense region of the 
prime halo. 
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Fig. 6. Snapshots of disk evolution and its stellar population. Stellar ages are 0-0.2 Gyr (blue), 0.2-2 Gyr 
(blueish), 2-5 Gyr (yellow) and older than 5 Gyr (red). In this model, a long-lived stellar bar is triggered 
by the interaction with an inclined prograde subhalo (hosting another disk) at z ~ 2.56 (left). Numerous 
interactions with additional subhalos (e.g., middle frames) leave the bar intact. A nuclear bar developes 
after z ~ 0.1, orthogonally to the prime bar (right). From Romano-Diaz et al. (2008c). 



Our main interest in subhalo populations 
lies in its impact on the disk evolution. 
Specifically, we focus on the comparison be- 
tween the PDM and BDM models. The envi- 
ronmental effects are clearly seen in the over- 
all decrease in A^h, both with and without 
baryons. This trend accelerates with z and es- 
pecially for the innermost massive PDM sub- 
halos. Hence we expect an excess of inner 
BDM subhalos over the PDM ones at lower 
redshifts. Fig. 5 shows that indeed the BDM 
subhalos dominate that inner 30 kpc region of 
the prime halo at low z, well after the epoch of 
major mergers which ends by z ~ 1.5 in these 
models. 

The time variable influx of the subhalos is 
clearly visible in Fig. 5. The maxima at z ~ 
4-2 (left frame) are associated with the ma- 
jor mergers. The subhalo population at z S 1 
(right frame) is much more steady, still one can 
observe 'waves' of subhalos that penetrate the 
inner 30 kpc. A more careful analysis shows 
that these subhalos have been accreted via fil- 
aments). They appear to cluster within the fil- 
ament but their merging there is suppressed. 
The ability of these subhalos to penetrate the 
central region of the prime halo is enhanced 
because they are 'glued' by the baryons and 
hence can withstand larger tides. 

Interactions between the DM substructure 
and stellar bars can be classified very roughly 
into flybys and direct hits (those include merg- 
ers). Prograde encounters can excite bars, ret- 
rograde can only cause lopsidedness, while di- 



rect hits weaken existing bars. Like the canon- 
ical bars formed as a result of the intrinsic 
instability, tidally-induced bars are very re- 
silient. They avoid a number of difficulties as- 
sociated with the linear stage of the bar in- 
stability. An example of such a long-lived bar 
triggered by a DM subhalo on a prograde in- 
clined orbit is shown in Fig. 6. Interactions 
with the surrounding substructure leave it in- 
tact (Gauthier et al. 2006; Romano-Diaz et 
al. 2008c). Without much exaggeration, one 
can say that in order to destroy this bar, one 
should largely destroy the disk which hosts it, 
although it does go through phases when it is 
quite weak. 

To summarize, there is a clear advantage 
to processes which involve tidal triggering of 
galactic bars, at least at high to intermediate, 
and possibly low z. This triggering can result 
either from an asymmetric DM distribution or 
from interactions with the surrounding sub- 
structure. Such a triggering constitute a finite 
amplitude perturbation and is independent of 
a number problems associated with the canon- 
ical (spontaneous) bar instability. Recent nu- 
merical simulations involving baryons show a 
more dominant population of subhalos in the 
innermost DM halo, in the regions which host 
growing galactic disks, and after the epoch 
of major mergers. Finally, the availability of 
substructure around field and cluster galaxies 
points to the bar fractions being (reasonably) 
independent of the environment. 
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